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respiration  rates  and  substrate  types. 

ABSTRACT 

Bacteria  of  the  group  “Dehalococcoides”  display  the  ability  to  respire  recalcitrant  chlorinated 
organic  compounds.  “Dehalococcoides”  strains?  respiratory  pathways  and  function  of  most 
genome-encoded  enzymes  responsible  for  dechlorination,  reductive  dehalogenases  (RDases), 
remain  incompletely  annotated.  To  further  the  description  of  the  biological  organization  of 
“Dehalococcoides”,  this  study  monitored  the  trancriptomic  response  of  “Dehalococcoides 
ethenogenes”  stain  195  using  two-color  microarrays.  This  study  analyzed  the  transcriptome  of  23 
varied  continuous  feed  (or  pseudo-steady  state  (PSS))  conditions  and  two  distinct  batch  fed 
conditions.  The  continuous  feed  experiments  were  comprised  of  57  cultures  with  varying  electron 
acceptor  feed  rates  (0-504  ?eeq/(L-hr)),  electron  acceptor  types  (tetrachloroethene  (PCE), 
trichloroethene  (TCE),  dichloroethene  (DCE),  2,3-dichlorophenol  (DCP),  and  no  electron 
acceptor),  electron  donor  to  acceptor  ratios  (0.7  to  17  on  an  electron  equivalence  (eeq)  basis),  and 
electron  donor  type  (butyrate,  lactate,  yeast  extract,  fermented  yeast,  pure  hydrogen,  or 
endogenous  biomass  decay).  When  similarly  respiring  (-120  ?eeq  PCE/(L-hr))  batch  and  PSS 
cultures  were  contrasted,  the  RDases  DET1545  and  DET0180  were  up-regulated  in  the  PSS 
cultures  indicating  activity  at  lower  overall  electron  acceptor  concentration.  For  all  continuous  fed  chloroethene 
cultures,  members  of  the  RDase  family  and  electron  transport  chain  displayed 
unique  clusters  of  transcripts  responding  either  positively,  negatively,  or  indifferently  to  the 
respiration  rate.  An  RDase  within  the  indifferent  group,  DET1 171,  was  highly  (3 1  ±  15fold)  up 
regulated  in  response  to  DCP  being  fed  as  the  electron  acceptor.  DET1 171  could  potentially  play 
a  role  with  DET1559  (27±4. 1  fold  up-regulated  during  DCP  growth)  and  DET03 1 8  pceA  in  the 
dechlorination  of  chlorophenols. 
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1  Abstract: 

2  Bacteria  of  the  group  “ Dehcilococcoides ”  display  the  ability  to  respire  recalcitrant  chlorinated 

3  organic  compounds.  “ Dehcilococcoides ”  strains’  respiratory  pathways  and  function  of  most 

4  genome-encoded  enzymes  responsible  for  dechlorination,  reductive  dehalogenases  (RDases), 

5  remain  incompletely  annotated.  To  further  the  description  of  the  biological  organization  of 

6  “Dehcilococcoides” ,  this  study  monitored  the  trancriptomic  response  of  “Dehcilococcoides 

7  ethenogenes ”  stain  195  using  two-color  microarrays.  This  study  analyzed  the  transcriptome  of  23 

8  varied  continuous  feed  (or  pseudo-steady  state  (PSS))  conditions  and  two  distinct  batch  fed 

9  conditions.  The  continuous  feed  experiments  were  comprised  of  57  cultures  with  varying  electron 

10  acceptor  feed  rates  (0-504  peeq/(L-hr)),  electron  acceptor  types  (tetrachloroethene  (PCE), 

11  trichloroethene  (TCE),  dichloroethene  (DCE),  2,3-dichlorophenol  (DCP),  and  no  electron 

12  acceptor),  electron  donor  to  acceptor  ratios  (0.7  to  17  on  an  electron  equivalence  (eeq)  basis),  and 

13  electron  donor  type  (butyrate,  lactate,  yeast  extract,  fermented  yeast,  pure  hydrogen,  or 

14  endogenous  biomass  decay).  When  similarly  respiring  (-120  peeq  PCE/(L-hr))  batch  and  PSS 

15  cultures  were  contrasted,  the  RDases  DET1545  and  DET0180  were  up-regulated  in  the  PSS 

16  cultures  indicating  activity  at  lower  overall  electron  acceptor  concentration.  For  all  continuous¬ 
ly  fed  chloroethene  cultures,  members  of  the  RDase  family  and  electron  transport  chain  displayed 

18  unique  clusters  of  transcripts  responding  either  positively,  negatively,  or  indifferently  to  the 

19  respiration  rate.  An  RDase  within  the  indifferent  group,  DET1 171,  was  highly  (31  ±  15fold)  up- 

20  regulated  in  response  to  DCP  being  fed  as  the  election  acceptor.  DET1 171  could  potentially  play 

21  a  role  with  DET1559  (27+4.1  fold  up-regulated  during  DCP  growth)  and  DET0318  pceA  in  the 

22  dechlorination  of  chlorophenols. 

23 

24 

25 

26 
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1  Introduction: 

2  Recalcitrant  halogenated  organics  contaminate  a  majority  of  active  hazardous  waste  cleanup  sites 

3  (55).  Site  managers  and  engineers  have  turned  toward  utilizing  populations  of  the  bacteria 

4  ‘ ‘ Deh alococco ides  ”  (DHC)  in  bioremediation  efforts  (47).  DHC  strains  are  capable  of 

5  anaerobicallyreductively  dehalogenating  halogenated  compounds  including  chlorinated-ethenes 

6  (48.50.54),  -ethanes  (18,25,54).  -phenols  (2,21),  -benzenes  (1,3,21),  -dioxins  (11,13,22),  - 

7  dibenzodioxins  (21),  -dibenzofurans  (46),  -naphthalene  (21),  and  -biphenyls  (8,21,80)  allowing 

8  for  flexibility  in  the  bioremediation  application  of  this  bacteria.  Additionally,  DHC  is  often  found 

9  and  active  at  chlorinated  organic  contaminated  sites  (7,8,10,12,14,24,31,35,44,81-82)  and  need 

10  only  to  be  biostimulated  to  begin  remediating  the  chlorinated  organic  compounds.  Several  case 

11  studies  focusing  on  bioaugmentation  of  DHC  at  PCE  and/or  TCE  contaminated  sites  have  also 

12  showed  success  (19,45,52,67). 

13  The  first  strain  of  DHC  identified,  “ Dehalococcoides  ethenogenes  strain  195”  (DET),  was 

14  isolated  in  1997  (53)  and  its  genome  sequenced,  assembled,  and  annotated  in  2005  (69).  Genomes 

15  of  additional  DHC  members,  such  as  CBDB1  (42,74),  VS  (58,75,49),  GT  (71),  and  BAV1 

16  (27,40,49),  have  been  sequenced  and  annotated  since,  with  other  members  being  isolated  and/or 

17  characterized  (e.g.  strain  FL2(28),  KB1(17)  metagenome  available  through  the  Joint  Genome 

18  Institute).  However,  much  of  the  biological  organization  of  DHC  remains  unknown.  DHC  is  an 

19  obligate  hydrogenotroph  and  dehalorespirer  (32).  The  identified  (69)  putative  enzymes  and 

20  regulatory  elements  comprising  the  electron  transport  chain  include  a  unique  ensemble  of 

21  reductive  dehalogenases  (RDases),  which  reduce  halogenated  organics;  hydrogenases,  which 

22  oxidize  DHC’s  sole  electron  donor  (hydrogen);  and  other  oxidoreductases  of  unknown  or 

23  misannotated  function  (e.g.  Fdh,  format  dehydrogenase;  Nuo,  NADH  ubiquinone  oxidoreductase; 

24  Mod,  molybdopterin  containing  oxidoreductase).  Every  DHC  strain  has  an  unique  subset  of  over 

25  90  RDases  (41)  likely  imparting  capabilities  for  dehalogenating  different  chlorinated  substrates  in 

26  unique  patterns  (33,34,51,79).  Of  DET’s  seventeen  functional  RDases,  DET0079  (TceA), 
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1  DET0318  (PceA),  DET1559,  and  DET1545  are  the  most  investigated  across  multiple 

2  experiments  within  a  mixed  community  culture  containing  DET  as  the  sole  dehalogenator  (Donna 

3  II)  (23,36,37,43,51,60-62,75).  DET’s  NiFe  hydrogenase  Hup  is  the  only  hydrogenase  predicted  to 

4  face  the  periplasm  (56,69,73).  Although  RDases,  Hup,  and  other  membrane  bound 

5  oxidoreductase  proteins  have  been  documented  in  several  DHC  studies  (57,77),  the  complete 

6  organization  and  interaction  between  these  respiration  elements  is  not  yet  fully  understood  (59). 

7  Observations  of  broad  expression  patterns  across  many  experiments  generate  hypothesizes  about 

8  which  transcripts  are  working  in  concert  to  fulfill  specific  biological  functions. 

9  To  broaden  the  understanding  of  the  biological  organization  of  “Dehcilococcoides”,  this  study 

10  employed  microarray  technology  to  monitor  the  transcriptomic  response  of  “Dehalococcoides 

11  ethenogenes  stain  195”  in  a  mixed  community  culture  to  19  varied  experimental  continuous-feed 

12  or  pseudo  steady  state  (PSS)  conditions  (a  total  of  57  total  cultures  are  monitored  considering 

13  biological  replicates).  In  these  studies,  PSS  conditions  refer  to  a  steady  mRNA  transcript  level 

14  and  not  biomass.  Previous  transcriptomic  studies  performed  on  DHC  include  characterizing 

15  unknown  species  or  cultures  (16,78)  and  monitoring  both  pure  and  mixed  well  defined  cultures 

16  (37,76).  The  microarray  studies  investigating  the  entire  transcriptome  under  pure-culture  and 

17  batch  fed  conditions  with  trichloroethene  (TCE)  and  hydrogen  fed  (as  the  electron  acceptor  and 

18  electron  donor  respectively)  (37)  revealed  key  insights  into  the  organism’s  time-dependent 

19  response,  carbon  metabolism  (72),  corrinoid  response  (38),  and  predicted  nitrogen  fixation 

20  pathway  (43). 

21  In  order  to  extend  the  existing  microarray  datasets  and  more  closely  simulate  field  conditions  this 

22  study  analyzed  the  transcriptome  of  cultures  grown  in  conditions  with  varying  electron  acceptor 

23  feed  rates  (from  0-504  peeq/(L-hr)),  electron  acceptor  types  (between  tetrachloroethene  (PCE), 

24  trichloroethene  (TCE),  dichloroethene  (DCE),  2,3-dichlorophenol  (DCP),  and  no  electron 

25  acceptor),  electron  donor  to  acceptor  ratios  (from  0.7  to  17  on  an  electron  equivalence  basis),  and 

26  electron  donor  type  (between  butyrate,  lactate,  yeast  extract,  fermented  yeast,  pure  hydrogen,  or 
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1  endogenous  biomass  decay  only).  By  continuously  feeding  substrates  over  the  course  of  days,  a 

2  near  constant  (PSS)  pool  of  mRNA  was  established.  A  large  enough  dataset  potentially  provides 

3  the  association  of  specific  genes’  perturbations  to  gradations  in  experimental  conditions  rather 

4  than  responses  confounded  by  non-steady  conditions  or  batch  growth.  By  monitoring  the 

5  transcriptome,  the  genome- wide  response  can  be  described.  Additionally,  the  growth  of  DET  in  a 

6  non-pure  culture  setting  allows  a  fuller  insight  into  the  maintained  mixed  community  cultures 

7  (Cornell  Donna  II  (53),  Pinellas  (26),  Victoria  (30),  KB-1  (18)  and  ANAS  cultures  (64))  and  a 

8  better  representation  of  field  conditions  in  which  DET  will  compete  with  methanogens  for 

9  hydrogen  produced  from  organisms  that  ferment  organic  substrates,  when  provided  (20,65,70). 

10  Materials  and  Methods: 

11  Culture  Growth  Conditions 

12  All  subcultures  studied  in  this  experiment  were  100  mL  samples  from  a  maintained  6L  mixed 

13  community  culture  (Donna  II)  previously  described  (20,60-62,65)  The  culture  contains  a  DET 

14  population  (near  5xl0x  cells/mL;  60%  of  total  population  on  per  cell  basis)  mixed  with 

15  methanogens  and  organisms  that  ferment  organic  substrates. 

16  Continuous  Feed  or  Pseudo-Steady  State  Setup  (PSS): 

17  The  continuous  feed  or  pseudo-steady  state  (PSS)  culture  setup  has  been  previously  described 

18  (60).  In  brief,  160  mL  glass  serum  bottles  were  filled  with  100  mL  of  3 -day  starved  mixed 

19  community  culture  sampled  from  the  main  reactor  (Donna  II).  Varying  combination  of  Pressure 

20  Lok  syringe  volumes  (10  mL,  5  mL,  2.5  mL,  500  pL,  and  100  pL;  VICI  Precision  Sampling)  and 

21  concentrations  of  electron  acceptor  and  donor  in  media  were  loaded  on  a  Cole  Palmer  74900 

22  Syringe  Pump  to  deliver  an  expected  continuous  rate. 

23  Gas  Chromatography  Monitoring 

24  Chlorinated-ethenes,  ethene,  and  methane  were  monitored  utilizing  a  gas  chromatography  method 

25  previously  described  (62).  Chlorinated  phenols  were  monitored  utilizing  a  modified  method 

26  previously  outlined  (39).  The  chlorinated  phenol  samples  were  comprised  of  1  mL  liquid  samples 


5 


1  withdrawn  from  the  subcultures  and  frozen.  The  thawed  samples  were  transferred  to  a  1  mL  glass 

2  vial,  100  uL  of  hexane  was  added  to  the  sample,  and  the  vial  was  shaken  at  300  lpm  for  six  hours. 

3  The  FID  gas  chromatography  method  selected  has  an  injection  temperature  of  200°  C,  a  constant 

4  column  temperature  of  160°  C,  and  the  FID  detector  temperature  of  200°  C  on  a  Supelcowax  10 

5  column  (30  m  x  0.53  mm  with  1.00  pm  film,  Supelco).  The  gas  mixture  was  200  kPa  nitrogen  as 

6  the  carrying  gas,  200  kPa  air,  and  250  kPa  hydrogen.  5  pL  of  the  sample  in  hexane  was  injected. 

7  Chemical  standard  curves  were  made  in  both  hexane  and  media  from  pure  chlorinated  phenols 

8  (ACROS  Organics). 

9  Calculation  of  respiration  rates 

10  Respiration  is  considered  to  occur  for  the  higher  chlorinated  ethenes  to  VC  or  ethene  (60)  with  6, 

11  4,  and  2  electron  equivalents  (eeqs)  for  PCE,  TCE,  and  DCE  respectively  with  the  specific 

12  calculation  previously  published  (61).  Extending  this  analysis  to  the  dichlorophenol  DCP 

13  experiments,  monochlorophenol  (MCP)  is  the  terminal  product  (21).  Therefore,  DET  will  receive 

14  2  eeqs  per  mol  DCP.  The  respiration  rate  equation  for  DCP,  rDCp,  in  peeq/L/hr  is: 

d  (MCP) 

r°cp  ~  2^r~ 

15  Microarray  Design: 

16  The  microarray  designed  for  this  experiment  was  an  Agilent  Technologies©  two-color,  15k,  60 

17  mer,  8  plex  array.  The  specific  designs  of  the  probes  utilized  a  modified  method  provided  by  the 

18  eArray  ©  software  suite  (5).  The  first  probe  set  utilized  for  the  batch  versus  PSS  experiment 

19  employed  a  base-composition  (BC)  technique  for  designing  and  scoring  the  best  probe  for  each 

20  transcript  (5).  The  second  probe  set  redesigned  for  the  remainder  of  the  experiments  was 

21  expanded  to  include  community  members,  non-protein  encoding  RNA  transcripts  (rRNAs, 

22  tRNAs),  and  a  luciferase  control.  A  modified  temperature  matching  (Tm)  method  developed  the 

23  probe  set  for  redesigned  array.  The  Tm  method  searched  for  an  optimal  design  with  a  consistent 

24  melting  point  temperature  (80°  C)  while  sacrificing  the  overall  quality  of  the  individual  probe  (5). 
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However,  if  a  BC  score  of  3,4,  or  poor  was  reported  for  a  transcript,  multiple  probes  around  the 
melting  temperature  of  80°  C  were  designed.  The  probe  with  the  best  base  composition  score 
nearest  to  the  80°  C  temperature  was  selected.  The  designs  were  searched  using  the  Basic  Local 
Alignment  Search  Tool  (BLAST)  (15)  against  both  the  National  Center  for  Biotechnology 
Information  (NCBI)  nucleotide  collection  and  the  assembled  mixed  community  metagenome 
(metagenome  data  not  reported,  currently  being  compiled  by  the  Joint  Genome  Institute)  to 
confirm  the  specificity  of  all  probe  sets.  Both  microarray  platform  designs  are  uploaded  and 
freely  available  at  the  NCBI  Gene  Expression  Omnibus  (GEO)  database 
(http://www.ncbi.nlm.nih.gov/geo/). 

RNA-cDNA  Handling  for  Microarray  Monitoring 

50  mL  of  liquid  culture  samples  were  centrifuged-at  14190xg.  The  centrifuged  sample  was  split 
into  8  individual  RNA  extractions  with  each  sample  following  the  RNeasy©  Mini  Kit  (Qiagen) 
extraction  previously  outlined  (62).  The  8  distinct  RNA  extractions  were  recombined  on  the  spin 
filter  before  the  first  RW1  buffer  wash.  The  Superscript  I  ©  DNAse  RNA  cleanup,  amino-allyl 
cDNA  formation,  cDNA  cleanup,  and  cDNA  labeling  with  Cy3  or  Cy5  followed  the  method 
outlined  (76).  The  quality  and  quantity  of  the  RNA  was  determined  using  the  RNA  6000  Nano 
assay  on  an  Agilent  2100  bioanalyzer  (Agilent  Technologies).  The  quantity  of  resulting  cDNA 
was  determined  by  using  the  Quant-IT™  OliGreen®  ssDNA  Assay  Kit  (Invitrogen).  A  common 
control  RNA  pool  sampled  from  the  main  Donna  II  reactor  after  3  days  of  starvation  was  labeled 
with  Cy3. 

Microarray  Hybridization  and  Scanning: 

For  each  experiment,  Cy5  labeled  cDNA  from  the  mixed  community  mRNA  pool  was  hybridized 
against  an  aliquot  of  common  control  of  Cy3  labeled  cDNA  from  3 -day  starved  culture.  The 
hybridization,  washing,  and  scanning  of  the  microarray  samples  was  performed  by  the  Cornell 
University  Microarray  Core  Facility  (http://cores.lifesciences.cornell.edu/brcinfo/)  and  followed 
the  methods  outlined  by  the  manufacturer  (5).  The  general  procedure  mixed  25  pi  (-400  ng)  of 
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1  the  labeled  cDNA  sample  with  25  pi  2x  Gene  Expression  (GEx)  Hybridization  Buffer  HI-RPM 

2  (5),  hybridized  the  sample  to  the  microarray  slide  at  65°  C  for  17  hours,  washed  with  GEx  Wash 

3  Buffer  1  and  2  (5)  at  room  and  elevated  (37°  C)  temperatures,  and  scanned  with  an  Agilent 

4  Technologies  Scanner  G2505C  with  a  5  pm  resolution. 

5  Statistical  Treatmen  t  of  the  Data  Set: 

6  Microarray  image  analysis  was  conducted  using  Agilent  Feature  Extraction  10.5  Image  Analysis 

7  Software.  The  Feature  Extraction  Software  was  also  utilized  to  perform  a  within  array  modified 

8  LOESS  normalization  between  the  Cy5  and  Cy3  signals,  to  calculate  a  log  ratio  between  the  Cy5 

9  and  Cy3  channels,  and  to  calculate  a  modified  Student  t-test  p-value  between  the  Cy5  and  Cy3 

10  signal  distributions  (83).  A  between  array  normalization  was  not  performed  due  to  the  wide 

11  variety  of  experimental  conditions  analyzed,  possibly  violating  the  assumption  that  less  than  20% 

12  of  transcripts  are  differentially  regulated  throughout  the  experimental  set.  The  more  detailed 

13  treatment  the  Agilent  Feature  Extraction  employed  can  be  found  in  the  user  manual  (4).  Replicate 

14  spots  for  the  same  probe  (ranging  from  6-20  spots/probe)  were  geometrically  averaged. 

15  Hierarchical  clustering  was  performed  in  R  2. 1 1 . 1  (http://www.r-project.org/).  The  K- Means 

16  cluster,  figure  of  merit  (FOM)  test,  and  the  Cluster  Affinity  Search  Technique  were  performed  in 

17  the  Multiple  Experiment  Viewer  (MeV  4.6)  freely  available  from  JCVI  (66).  The  raw  and 

18  normalized  data  is  freely  available  at  the  NCBI  GEO  database. 

19  Results: 

20  Batch  Growth  Compared  to  the  Continuous  Feed  Pseudo  Steady  State  System 

21  The  first  set  of  microarrays  compared  the  mixed  community  culture  with  similar  electron 

22  acceptor/donor  conditions  under  batch  and  pseudo-steady  state  growth.  The  batch  culture 

23  dechlorinated  660  peeqs/L  of  PCE  to  VC  or  ethene  in  6  hours  (110  peeqs/(L-hr))  while  the  PSS 

24  experiments  received  a  continuous  input  of  53+13  or  83+2.8  peeqs/L-hr  of  PCE  or  TCE 

25  respectively.  The  main  biological  differences  are  the  high  initial  pulse  of  electron  acceptor  and 

26  donor  seen  by  the  batch  cultures  to  the  low  concentrations  maintained  in  the  PSS-setup  (the  batch 
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1  culture  experiences  the  initial  pulse  of  660  peeqs/L  while  the  PSS  cultures  maintain  a  constant 

2  level  of  about  12+9.5  peeqs/L  of  chlorinated  ethenes  based  on  the  last  48  hours  of  GC  data)  and 

3  that  the  rate  of  continuous  feeding  in  the  PSS  cultures  was  maintained  for  7  days.  The  two  PCE  (1 

4  &  2)  and  two  TCE  (1  &  2)  experiments  are  compared  via  a  correlogram  (Figure  1).  The  PCE  (1  & 

5  2)  cultures  do  not  display  a  high  correlation  with  each  other  due  to  the  inadvertent  variation  in  the 

6  feeding  of  electron  acceptor  (PCE  land  2  were  fed  66  and  44  peeq/(L-hr)  respectively,  (GC  data 

7  displayed  in  Supplemental  Image  1).  To  determine  the  transcripts  differentially  regulated 

8  similarly  across  all  four  arrays,  only  transcripts  where  all  replicate  probe  spots  (n=6-20)  on  the 

9  microarray  exceed  a  fold  ratio  between  the  PSS  and  batch  cultures  of  two  (for  up-regulated 

10  genes)  or  one-half  (for  down-regulated  genes)  and  are  below  a  p-value  of  0.05  are  considered. 

11  Overall,  compared  to  the  batch  control,  PCE  1,  PCE  2,  TCE  1,  and  TCE  2  displayed  193,  165, 

12  144,  and  275  differentially  regulated  transcripts,  respectively.  Of  these,  43  (23  up-  and  20  down- 

13  regulated)  are  shared  across  all  samples  and  displayed  in  Table  1.  Within  the  23  up-regulated 

14  transcripts,  the  RDase  DET0180  and  its  anchoring  protein  as  well  as  the  RDase  anchoring  protein 

15  DET1544  display  a  5. 5+4.1 1,  5.1+0.92,  and  3.5+0.59  fold  increase  respectively.  The  presence  of 

16  the  anchoring  protein  DET1544  without  the  catalytic  RDase  DET1545  transcript  is  suiprising  as 

17  they  are  predicted  to  reside  on  the  same  polycistronic  mRNA.  All  probes  for  the  RDase  DET1545 

18  pass  the  requirement  of  displaying  a  p-value  below  0.05,  but  the  transcript  is  excluded  from  the 

19  stringent  analysis  due  to  two  probe  spots  (out  of  20  replicates)  in  the  PCE1  experiment  failing  to 

20  exceed  the  ratio  cut  off  of  2.  If  the  stringency  of  the  analysis  is  relaxed  to  remove  the  log  ratio 

21  cutoff,  the  number  of  transcripts  shared  across  all  experiments  is  increased  from  43  to  65 

22  (Supplemental  Table  1).  DET1545  falls  within  this  category  displaying  an  average  ratio  across  all 

23  experiments  of  3.3+0.86.  The  20  significantly  down-regulated  transcripts  under  PSS  conditions  in 

24  the  stringent  analysis  are  populated  by  many  key  transcripts  encoding  for  enzymes  in  respiratory 

25  energy  conservation.  These  include  transcripts  of  the  ATPcise  operon  (DET0558,  DET0560, 

26  DET0561,  DET0562)  and  subunits  of  the  hydrogenase  ech  (DET0864).  Additional  members  of 
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1  the  ech  operon  (DET0861  and  DET0862)  are  found  in  the  expanded  data  set  (Supplemental  Table 

2  1). 

3  Expression  profiles  comparison  across  multiple  PSS  experiments: 

4  In  order  to  explore  a  wide  variety  of  environmental  conditions,  mRNA  pools  from  53 

5  experiments  (performed  on  14  separate  days)  were  analyzed  by  the  microarray  technique 

6  described  above.  The  parameters  and  conditions  for  each  experiment  are  outlined  in  Table  2.  The 

7  electron  acceptors  were  varied  from  none,  PCE,  TCE,  DCE,  and  DCP  while  the  electron  donor 

8  for  the  PCE  experiments  was  varied  from  none,  butyrate  (with  yeast  extract  (YE)),  lactate  (with 

9  YE),  hydrogen,  YE,  and  fermented  yeast  extract  (FYE).  The  feed  rates  of  electron  acceptor  and 

10  the  ratio  of  electron  acceptor  to  donor  were  altered  as  well.  The  respiration  rates  varied  between  0 

11  peeq/(L-hr)  in  cultures  not  fed  electron  acceptor  to  170  peeq/(L-hr)  in  the  PCE  High 

12  Experiments. 

13  When  the  experiments  are  organized  via  a  hierarchal  Ward  cluster  dendogram  based  on  a 

14  Euclidean  distance  matrix  of  the  log  ratio  fold  changes  (Figure  2),  two  distinct  groupings  develop. 

15  The  two  clusters  are  separated  primarily  by  the  respiration  rate  recorded  with  higher  and  lower 

16  respiration  rates  clustering  together.  However,  the  High  2  from  the  High-Low  PCE  experiment 

17  (HLH2_INHIB)  plots  with  the  higher  respiration  cluster  when  in  fact  respiration  ceased  prior  to 

18  culture  sampling  due  to  the  formation  of  non-aqueous  phase  PCE.  HLH1_INHIB  and 

19  HLH2_INHIB  display  a  unique  decay  series  as  both  cultures  experienced  an  initial  high 

20  respiration  rate  (-100  peeq/(L-hr) ),  but  PCE  built  up  above  its  solubility  limit  at  day  2  inhibiting 

21  PCE  respiration.  Additionally,  the  DCE  “3  rate”  mid-range  (DCE  B 1  and  B2)  do  not  fall  into  the 

22  same  cluster.  The  remainder  of  the  experiments  and  biological  replicates  plot  as  expected.  The 

23  hydrogen  abiotic  controls  (H2A1  and  A2)  cluster  with  the  3  day  decay  experiment  (DecayAl)  as 

24  both  are  decay  experiments  on  the  same  time  scale.  Additionally,  the  experiments  fed:  PCE  with 

25  either  no  electron  donor,  YE,  or  FYE  with  PCE  (P0FY_XX);  the  chlorophenol  experiment  (Cp); 
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1  the  PCE  low  series  (PLL);  the  decay  7-day  (DecayB)  replicates;  the  PCE  high  series  (HiP);  the 

2  half  butyrate  series  (PHB);  and  the  majority  of  the  “3  rate”  experiments  (P3,  T3,  and  D3)  cluster 

3  with  their  respective  conditions  and  biological  replicates,  supporting  the  biological 

4  reproducibility  of  the  overall  transcriptional  response. 

5  Gene  clusters 

6  After  examining  the  clustering  of  experiments,  the  clustering  of  genes  with  similar  behavior  was 

7  explored  with  a  K-means  clustering  technique.  This  technique  was  able  to  align  all  genes  into 

8  biologically  relevant  clusters  (63).  However,  the  hindrance  to  employing  a  K-means  clustering 

9  technique  was  that  the  number  of  clusters  assigned  must  be  set  prior  to  analysis.  Based  on  a  figure 

10  of  merit  test,  the  number  of  biological  clusters  was  set  at  20.  A  hierarchal  dendrogram  displaying 

11  the  relationship  among  the  twenty  clusters  is  displayed  in  Figure  3  with  the  specific  gene  list  of 

12  each  cluster  described  in  Supplemental  Table  2.  Other  clustering  techniques,  such  as  the  cluster 

13  affinity  search  technique  (CAST)  (9),  gave  similar  results  (data  not  shown).  The  20  clusters  map 

14  into  four  distinct  superclusters.  Figure  3.b-e  displays  the  cluster’s  centroid  log  ratio  expression 

15  value  versus  respiration  rate  broken  into  the  four  primary  clusters.  Figure  3.b  displays  a  general 

16  positive  trend  in  response  to  respiration  rate  while  Figure  3.c  displays  a  general  negative  trend 

17  (Figure  3.d-e  displays  clusters  with  centroids  of  these  clusters  remain  nearly  constant).  This 

18  bifurcation  is  important  to  note  as  the  majority  of  high  respiration  rate  linked  RDases  ( tee  A ,  pceA, 

19  DET1559)  and  respiration  chain  elements  {hup,  ATPcise )  are  members  of  the  clusters  plotted  in 

20  Figure  3.b  while  the  low  respiration  linked  RDases  (DET0173,  DET1535,  DET1538,  and 

21  DET1545)  are  members  of  the  clusters  plotted  in  Figure  3.c. 

22  Across  all  clusters,  three  centroid  data  points  occurring  around  the  35  peeq/(L-hr)  respiration  rate 

23  appear  to  be  outliers  (lower  in  Figure  3.b  and  higher  in  Figure  3.c).  These  centroids  are 

24  noteworthy  as  they  represent  the  DCE_INHIB  experiments.  If  the  final  48  hours  of  respiration 

25  prior  to  sampling  are  considered,  the  respiration  rate  was  actually  around  2.4±1 .8  peeq/(L-hr) 
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1  which  shifts  all  three  centroids  to  the  left  in  the  respiration  diagrams  of  Figure  3.b-e.  Following 

2  this  adjustment,  these  points  do  not  appear  to  be  outliers  but  support  the  overall  trend. 

3  K-Means  Clusters  Mapped  onto  DET’s  Genome: 

4  -Describe  binning  the  genes  based  on  supercluster 

5  -Show  that  there  are  several  regions  where  a  siginificant  number  of  genes  appear  in  the  same 

6  supercluster 

7  -Describe  the  virus  region  and  explore  other  genes  that  are  clustered  with  it. 

8  Average  Expression  Values  Across  all  Experiments 

9  -Describes  those  genes  that  are  most  highly  expressed 

10  -Compare  to  Johnson  et.  al  to  see  if  any  uniquely  turned  on 

11  Putative  Respiration  Chain  Element  Rate  Response 

12  The  fact  respiration  is  a  key  and  informative  phenotype  of  DET  prompted  interest  in  plotting 

13  expression  patterns  of  transcripts  versus  experimental  respiration  rate.  Focusing  on  the  responses 

14  of  mRNAs  encoding  for  known  or  predicted  respiration  chain  elements  could  determine  if  there 

15  are  unique  electron  transport  chain  members  responding  to  low  or  high  respiration.  Figure  4 

16  displays  the  expression  pattern  clustering  of  18  RDases  (17  full  length  plus  DET0162,  an  RDase 

17  with  a  known  point  mutation  (69)),  5  hydrogenases  ( hym ,  hyc,  ech,  hup,  vhu ),  other  membrane- 

18  bound  oxidoreductases  ( mod  (potential  molybdopterin  oxidoreductases ),fdh  (annotated  formate 

19  dehydrogenase),  nuo  (NADH-ubiquinone-oxidoreductase)),  and  ATPase.  In  this  figure  a  single 

20  polypeptide  subunit  is  selected  to  represent  the  full  operon.  However,  members  of  the  electron 

21  transport  chain  often  are  encoded  on  multi-subunit  operons.  A  version  of  Figure  4  with  all 

22  subunits  plotted  is  presented  in  Supplemental  Figure  2.  The  majority  of  subunits  of  the  respective 

23  operons  cluster  closely  together.  However,  variability  is  seen  in  the  nuo,  ech,  and  hym  transcripts. 

24  In  Figure  4,  three  groups  emerge  and  again  vary  in  their  trend  relative  to  respiration  rate.  The  first 

25  grouping  contains  the  high-respiration-linked  RDases  with  the  majority  of  the  other  predicted 

26  respiratory  elements.  As  seen  by  Figure  4.b,  these  transcripts  respond  very  sharply  to  increased 
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1  levels  of  respiration  and  then  hit  a  plateau.  This  plateau  was  noted  in  previously  published  work 

2  by  Rahm  and  Richardson  (60.61).  The  low-respiration-linked  RDases,  DET0173,  DET1535, 

3  DET1538,  and  DET1545,  are  tightly  clustered  in  the  second  group  and  at  rates  above  5  peeq/(L- 

4  hr)  display  an  inverse  relationship  to  respiration  as  first  noted  in  the  K-means  analysis  (Cluster  1 

5  in  Figure  3.c).  The  remaining  RDases  and  the  putative  respiration  chain  elements  mod  and  ech 

6  appear  to  fall  within  a  cluster  that  does  not  respond  greatly  to  varied  levels  of  chlorinated  ethene 

7  respiration  although  clustering  between  modA  (DET0103)  and  DET0180  is  observed  (Figure  4.a). 

8  Response  of  Respiration  Elements  to  Switching  the  Electron  Acceptor  Type  from  Chlorinated 

9  Ethenes  to  2,3 -Diehl orophenol  and  1,2, 3,4  Tetrachlorobenzene 

10  The  electron  acceptor  was  varied  in  the  experimental  set  among  PCE,  TCE,  DCE  or  2,3- 

11  dichlorophenol  (DCP).  Figure  5. a  displays  the  natural  value  differential  expression  of  the 

12  individual  RDases  across  a  relatively  consistent  low-rate  of  respiration  with  varying  the  electron 

13  acceptor  type  (electron  donor  was  always  butyrate).  The  tceA  transcript  is  about  equally  up- 

14  regulated  in  the  DCP  samples  as  with  any  of  the  chlorinated  ethenes  (PCE,  TCE,  DCE)  (Figure 

15  4.a).  Additionally,  DET1 171  and  DET1559  are  both  highly  up-regulated  in  the  chlorophenols  fed 

16  cultures  (32  ±16  and  27  +  4.1  fold  respectively).  The  distribution  of  the  absolute  intensity  of 

17  RDase  spots  for  the  varied  electron  acceptor  types  (Figure  5.b)  reveals  that  for  the  TCE  and  DCE 

18  fed  cultures  tceA  dominates  the  percentage  of  the  RDase  transcriptional  pool,  for  the  PCE  fed 

19  culture  DET0162  dominates,  while  for  the  DCP  fed  culture  pceA  is  the  dominant  transcript 

20  comprising  above  40%  of  the  total  RDase  processed  fluorescence.  Additionally,  pceA’s 

21  differential  expression  displays  a  slight  non-significant  up-regulation  for  the  DCP  culture  as 

22  compared  to  the  chlorinated  ethene  fed  cultures.  Furthermore,  DET0180  sees  a  decrease  in  the 

23  percentage  of  RDase  intensity  between  DCP  and  TCE/DCE. 

24  Including  the  other  respiratory  chain  elements  in  the  analysis  (data  not  shown)  reveals  that  the 

25  majority  of  the  transcripts  are  down  regulated  under  DCP  conditions.  However,  nuo  (DET1571), 
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1  vim  (DET0615),  and  ech  (DET0860)  do  not  display  a  significant  difference  between  the  DCP  fed 

2  and  chlorinated  ethene  fed  cultures  in  terms  of  ratio  to  standard  controls. 

3  The  analysis  for  DCP  specific  response  was  expanded  to  include  all  transcripts  in  a  rank  ordered 

4  list  based  off  of  differential  expression  ratios  for  the  DCP  experimental  series  (Supplemental 

5  Table  3).  The  first  entry  is  DET1 171,  as  previously  indicated  in  Figure  5. a.  Besides  the  two 

6  RDases  DET1 171  and  DET1559  reside  in  the  top  twenty  of  all  genes  by  fold  change,  response 

7  regulators  and  hypothetical  proteins  heavily  populate  the  list.  Two  histidine  kinase  sensory  boxes 

8  (DET0301  and  DET1560),  pyrR  (DET1 198,  putative  transcriptional  binding  protein  of 

9  pyrimidine  nucleotide  biosynthesis),  and  recX  (DET1607,  putative  regulatory  protein  for  the 

10  DNA  recombinase  RecA)  all  display  an  up  regulation  greater  than  ten-fold.  Among  the  most 

11  down-regulated  transcripts,  an  over  abundance  of  genes  related  to  iron  or  cobalamin  are  present. 

12  These  genes  include  two  ABC  transporters  (DET0650  (corrinoid  specific)  and  DET1 176)  and  a 

13  biotin-Acetyl-CoA-carboxylase  ligase  (DET0849,  also  putatively  described  as  a  transcriptional 

14  repressor). 

15  Tetrachlorobenzene  experiments  (two  weeks) 

16  Tetrachlorobenzene  arrays  (one  week) 

17  Analysis  (few  days) 

18  Discussion: 

19  In  the  comparison  of  the  batch  fed  culture  to  the  PSS  setup,  the  RDases  DET1545  and  DET0180 

20  are  displayed  to  be  highly  up-regulated  under  PSS  conditions.  This  suggests  that  DET1545  and 

21  DET0180  are  up-regulated  in  response  to  lower  electron  acceptor  availability  conditions 

22  (specifically  the  batch  fed  culture  processed  a  660  peeq/(L)  pulse  of  PCE  at  a  rate  of  1 10  peeq/(L- 

23  hr)  while  the  PSS  received  a  continuous  low  rate  of  55±8.2  peeq/(L-hr)),  as  supported  by 

24  previous  publications  (60-62).  Although  neither  transcript  is  dominant  in  the  absolute  intensity 

25  sense  (data  not  shown),  the  minor  RDases  potentially  serve  as  key  enzymes  when  electron 

26  acceptor  availability  is  low  or,  their  regulatory  elements  are  sensitive  to  external  electron  acceptor 
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1  concentrations.  When  the  analysis  is  expanded  to  include  putative  members  of  the  entire 

2  respiration  chain,  it  is  noted  that  the  ech  and  ATPcise  operons  are  both  down-regulated.  The 

3  ATPase  enzyme  encodes  for  a  proton  pump  utilized  primarily  for  energy  generation  and  the  Ech 

4  hydrogenase  putatively  encodes  for  a  proton  pump  which  generates  low  potential  electrons  that 

5  could  be  used  to  reduce  ferredoxin  or  possibly  cobalt  in  RDase  corrinoids  to  the  Co  I  state,  or  for 

6  other  biosynthetic  metabolites  (69).  Potentially,  the  batch  cultures  are  limited  by  the  rate  of 

7  hydrogen  production  from  butyrate  while  the  continuous-fed  cultures  are  limited  for  electron 

8  acceptor  therefore  allowing  higher  hydrogen  levels.  The  higher  hydrogen  level  in  the  continuous 

9  feed  cultures  may  mitigate  “ Dehcilococcoides ”  need  for  reverse  electron  transport  by  Ech. 

10  In  previous  work  with  the  mixed  community  culture  Donna  II,  qRTPCR  was  used  to  monitor 

11  expression  levels  of  selected  genes  (60-62).  The  trends  from  that  work  with  respect  to  the 

12  RDases’  relationship  to  respiration  rate  are  supported  here  with  genome  wide  profiles.  Namely, 

13  some  highly  expressed  respiration  chain  elements  and  RDases  up  regulate  sharply  in  response  to 

14  increased  respiration  rate  or  electron  acceptor  concentration  before  plateauing  at  the  highest  rates 

15  while,  in  contrast,  other  RDases  are  up-regulated  only  at  low  respiration  rates  or  electron  acceptor 

16  concentrations  before  being  down-regulated  at  higher  rates.  This  trend  is  also  supported  by  batch 

17  studies  with  DET  where  the  low-respiration  rate  linked  cluster  of  RDases  is  up-regulated  upon 

18  entry  into  stationary  phase  (37).  This  bifurcation  of  behavior  in  the  RDases  may  indicate  a  unique 

19  response  of  specializing  RDases  in  the  electron  transport  chain  and  prove  useful  for  the 

20  identification  of  field  site  biomarkers  in  determining  the  respiration  state  of  the  organisms  in  situ. 

21  However,  though  the  expression  patterns  of  these  RDases  are  remarkably  repeatable,  their 

22  substrate  range  remains  unknown  and  functional  studies  are  needed. 

23  Focusing  on  responses  specific  to  electron  acceptor  type,  previous  studies  comparing  PCE,  TCE, 

24  and  DCE  at  similar  respiration  rates  showed  shifts  in  relative  expression  of  the  RDases  being 

25  tracked  by  qRTPCR  (e.g.  pceA  was  down-regulated  slightly  when  TCE  or  DCE  was  fed),  but 

26  these  shifts  were  less  than  two  fold  (61).  Shifting  to  a  different  class  of  chloroorganic  election 
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1  acceptor  altogether  displayed  larger  shifts  (2,21,23).  The  current  study,  in  varying  the  electron 

2  acceptor  type  from  chlorinated  ethenes  to  DCP,  supports  the  majority  of  reported  findings  (23). 

3  DET1559  and  pee  A  have  been  highlighted  as  two  transcripts  potentially  involved  in  chlorophenol 

4  dechlorination  (23).  These  two  findings  are  supported  in  the  current  study  by  the  dominance  of 

5  pee  A  in  the  RDase  absolute  intensity  distribution  (Figure  5.b),  and  the  high  up-regulation  of 

6  DET1559  (27  ±4.1  fold  for  DCP  compared  to  4.4  ±1.4  fold  for  PCE).  PceA  has  been  shown  to 

7  be  bifunctional,  able  to  dechlorinate  both  chlorinated  ethenes  and  phenols  (2). 

8  Expanding  off  of  these  findings,  the  current  study  displays  that  the  RDase  DET1 171  is  highly  up 

9  regulated  in  the  DCP  cultures  (31  ±  15  fold)  while  not  significantly  up-regulated  in  the 

10  chlorinated  ethene  fed  cultures  (Figure  5. a).  In  the  Fung  et  al.  2007  study,  the  RDase  DET1 171  is 

11  indicated  as,  but  not  highlighted  for  due  to  low  overall  expression,  displaying  significant  up 

12  regulation.  When  DET1 171  ’s  amino-acid  sequence  was  searched  by  BLAST  against  the  NCB1 

13  protein  database,  ortholog  RDases  appear  in  the  CBDB1,  BAV1,  and  VS  strain  (95%,  95%,  and 

14  97%  amino  acid  identity  respectively)  as  previously  indicated  (49).  Additionally,  weaker  hits 

15  match  with  the  ortho-chlorophenol  RDases  in  the  Desulfitobacterium  genus  (39%  amino  acid 

16  identity)  with  seven  out  of  the  top  twenty  hits  described  as  a  chlorophenols-related  RDase. 

17  However,  the  closest  hits  outside  of  the  “ Dehalococcoides ”  group  are  genes  from  species  of 

18  Photobacterium  and  Shewanella,  which  may  be  non-energy  conserving  RDases  involved  in  the 

19  breakdown  of  a  chlororganics.  Therefore,  the  correlation  of  DET1 171  with  DCP  dechlorination 

20  may  be  the  result  of  shared  transcriptional  regulation  rather  than  actual  functional  dechlorination 

21  by  the  final  enzyme.  To  fully  explore  the  function  and  behavior  of  DET1 171  and  the  other  RDase 

22  family  members,  the  dataset  which  is  primarily  comprised  of  chlorinated  ethene  fed  cultures 

23  should  be  expanded  to  explore  more  perturbations  in  chlorophenols,  chlorobenzene, 

24  chlorodibenzodioxins,  and  other  chlorinated  substances. 
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9 

10 

11  Figure  1 .  Comparison  of  the  microarray  results  for  the  Batch  and  PSS  fed  cultures.  The 

12  correlogram  displays  the  similarity  of  experimental  replicates  TCE  1&2  and  PCE  1  &  2  by 

13  correlating  the  gene  expression  across  experiments  for  all  transcripts  compared  to  batch 

14  conditions. 

15 

16  Table  1.  The  43  genes  that  are  similarly  regulated  amongst  all  four  Batch  vs.  PSS  samples  with 

17  the  transcript  name,  regulation  (up  or  down  when  comparing  PSS  to  batch  level),  average  natural 

18  value  ratio  recorded,  and  the  standard  deviation  of  that  ratio  across  all  4  arrays  (40  spots  on  4 

19  arrays). 

20 

21  Table  2.  Experimental  conditions  for  the  continuous-feed  (PSS)  investigations.  The  experimental 

22  title  provides  a  brief  description  of  the  overall  experiment  where  the  culture  title  designates  the 

23  electron  acceptor  (PCE,  TCE,  DCE,  DCP,  or  none  (-) )  and  electron  donor  source  (butyrate, 

24  lactate,  hydrogen,  yeast  extract  (YE),  fermented  yeast  extract  (FYE),  hydrogen,  or  none  (-) )  type. 

25  Note,  the  experiments  that  received  butyrate  and  lactate  also  were  received  YE.  Additionally,  the 
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1  Batch  versus  PSS  direct  arrays  are  not  displayed  as  those  experiments  were  compared  against  a 

2  different  control. 

3 

4  Figure  2.  Full  experimental  Ward  cluster  dendogram  for  all  experiments  based  off  of  an 

5  Euclidean  distance  matrix  for  the  log  ratio  fold  change  of  all  gene  targets.  Consult  Table  2  for 

6  culture  descriptions. 

7 

8  Figure  3.  Breakdown  of  the  of  individual  genes’  expression  patterns  into  clusters,  (a)  A  Euclidean 

9  distance  and  Ward  cluster  dendogram  comparing  the  similarities  of  the  clusters,  (b-e)  Scatterplots 

10  plotting  the  centroid  value  (with  95%  confidence  intervals  about  the  mean)  against  the 

11  experimental  respiration  rates  for  all  chlorinated  ethenes. 

12 

13  Figure  4.  Respiratory  chain  elements’  response  to  respiration  rate,  (a)  Ward  hierarcial  clustering 

14  with  an  Euclidean  distance  matrix  of  the  differential  expression  in  log  ratio.  Transcripts  with  an 

15  (*)  do  not  exceed  the  intensity  cutoff  limit  of  100.  The  bootstrap  analysis  of  the  cluster  results  in 

16  two  scores,  the  approximately  unbiased  (AU)  p-value  in  red  and  the  bootstrap  probability  (BP)  in 

17  green,  (b-d)  Scatterplots  of  the  differential  expression  log  ratio  for  individual  genes  in  each  large 

18  cluster  plotted  against  the  respiration  rate  (peeq/L/hr).  The  black  point  represents  the  average  of 

19  the  log  ratio  values  in  each  gene  cluster  with  the  error  bars  representing  the  95%  confidence 

20  interval  about  the  mean. 

21 

22  Figure  5.  RDase  distribution  across  different  electron  acceptors  for  the  PCE,  TCE,  DCE,  DCP, 

23  and  No  Acceptor  fed  conditions  (a)  natural  value  ratio  compared  to  a  3-day  starved  culture  for  the 

24  17  RDase  and  the  rpoB  transcript  is  plotted.  Error  bars  represent  standard  error  across  biological 

25  replicates,  (b)  The  absolute  processed  intensity  value  distribution  of  the  RDases. 

26 
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Table  3.  Correlation  scores  of  RDases  that  meet  the  background  cutoff  to  neighboring  annotated 
genes  including  putative  response  regulators.  The  solid  black  line  separates  the  RDases  based  on 
their  clustering  and  response  to  respiration  rate  (RDases  responding  positively  to  respiration  rate 
above  the  black  line,  negatively  below)  displayed  in  Figure  4.  A  negative  correlation  score 
indicates  an  inverse  correlation.  Values  >0.25  or  <-.25  are  considered  significant  by  the  one- 
tailed  t-test  with  N=49. 


Table  4.  Correlation  scores  of  tceA  to  RDases  and  response  regulators  adjacent  to  RDases  (left 
column)  and  correlation  scores  of  tceA  to  all  transcripts  on  DET’s  genome  represented  on  the 


microarray  (right  column). 
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1  Figure  1.  Comparison  of  the  Batch  and  PSS  fed  cultures.  The  correlogram  displays  the  similarity 

2  of  experimental  replicates  TCE  1&2  and  PCE  1  &  2  by  correlating  the  gene  expression  across 

3  experiments  for  all  transcripts. 
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Figure  1.  Comparison  of  the  Batch  Control  and  PSS  fed  cultures  log  ratio  correlations.  The  correlogram 
displays  the  similarity  of  the  experimental  replicates  TCE  1&2  and  PCE  1&2  be  correlating  the  log  ratio 
for  every  gene  transcript  across  experiments. 
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1  Table  2.  Experimental  conditions  for  all  the  12  investigations  comparing  the  PSS  versus  a  3  day 

2  starved  control.  The  experimental  title  provides  a  brief  description  of  the  overall  experiment 

3  where  the  culture  title  designates  the  electron  acceptor  (PCE,  TCE,  DCE,  DCP,  or  ncne  (-)  )  and 

4  electron  donor  source  (butyrate,  lactate,  hydrogen,  yeast  extract  (YE),  fermented  yeast  extract 

5  (FYE),  hydrogen,  or  none  (-)  )  type.  Note,  the  experiments  that  received  butyrate  and  lactate  also 

6  were  treated  with  YE. 
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1  Figure  3.  Biological  breakdown  of  the  clusters,  (a)  A  Euclidean  distance  and  Ward  cluster 

2  dendogram  comparing  the  similarities  of  the  clusters,  (b-e)  Scatterplots  plotting  the  centroid  value 

3  with  95%  confidence  intervals  about  the  mean  against  the  experimental  respiration  rates  for  all 
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Predicted 

DNA 

Binding 

Response 

Reg. 

Predicted 

MarR 

Regulator 

Annotated 

Lipo¬ 

protein 

Other  Annotations 
(Including  Hypothetical) 

DET0079  -  tceA 

DET0078 

0.97 

DET0318 - pceA 

DET0319 

DET0315 

DET0316 

DET0321 

DET0317 

DET0320 

0.94 

0.18 

0.30 

-0.10 

-0.12 

0.69 

DET1559 

DET1558 

DET1560 

DET1561 

DET1557 

DET1556 

0.69 

0.60 

0.70 

0.61 

0.78 

DET0173 

DET0175 

DET0171 

DET0170 

DET0174 

DET0172 

DET0176 

0.95 

-0.25 

-0.19 

0.94 

0.39 

0.58 

DET0180 

DET0181 

DET0178 

DET0177 

DET0179 

0.68 

0.54 

0.48 

0.34 

DET1535 

DET1534 

DET1536 

DET1533 

DET1532 

0.94 

-0.63 

0.41 

0.49 

DET1538 

DET1537 

DET1539 

DET1540 

0.68 

0.48 

-0.10 

DET1545 

DET1544 

DET1541 

DET1542 

DET1543 

0.77 

0.33 

0.10 

-0.31 

1  Table  3.  Correlation  scores  of  RDases  that  meet  the  background  cutoff  to  neighboring  annotated 

2  genes  including  putative  response  regulators.  The  solid  black  line  separates  the  RDases  based  on 

3  their  clustering  and  response  to  respiration  rate  (RDases  responding  positively  to  respiration  rate 

4  above  the  black  line,  negatively  below)  displayed  in  Figure  4.  A  negative  correlation  score 

5  indicates  an  inverse  correlation.  Values  >0.25  or  <-.25  are  considered  significant  by  the  one- 

6  tailed  t-test  with  N=49. 
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1  Figure  4.  Respiration  chain  elements  group  response  to  respiration  rate,  (a)  Ward  hierarcial 

2  clustering  with  an  Euclidean  distance  matrix  of  the  differential  expression  in  log  ratio.  Transcripts 

3  with  an  (*)  do  not  exceed  the  intensity  cutoff  limit  of  100.  The  bootstrap  analysis  of  the  cluster 

4  results  in  two  scores,  the  approximately  unbiased  (AU)  p- value  in  red  and  the  bootstrap 

5  probability  (BP)  in  green,  (b-d)  Scatterplots  of  the  differential  expression  log  ratio  for  individual 

6  genes  in  each  large  cluster  plotted  against  the  respiration  rate  (ueeq/L/hr).  The  black  point 

7  represents  the  average  of  the  log  ratio  values  in  each  gene  cluster  with  the  error  bars  representing 

8  the  95%  confidence  interval  about  the  mean. 


mRNA  ID 

Correlation 

Score 

mRNA  ID 

Correlation 

Score 

DET0320  hypothetical  protein 

0.90 

DET1407  BNR/Asp  box  repeat  domain  protein 

0.98 

DET1559  RDase 

0.78 

DET0078  tceA  anchoring  protein 

0.97 

DET0170  DNA  binding  response  regulator 

0.76 

DET0511  hypothetical  protein 

0.93 

DET1543  hypothetical  protein 

0.69 

DET0112  [Ni/Fe]  hydrogenase 

0.93 

DET0171  sensory  box  sensor  histidine  kinase 

0.68 

DET1600  twin  arginine  translocation  TatB 

0.93 

DET0177  DNA  binding  response  regulator 

0.67 

DET0117  transcriptional  regulator 

0.92 

DET1556  hypothetical  protein 

0.67 

DET0352  hypothetical  protein 

0.91 

DET0319  RDase  anchoring  protein 

0.67 

DET0191  iron  sulfur  cluster  binding  protein 

0.91 

DET1557  lipoprotein,  putative 

0.63 

DET0320  hypothetical  protein 

0.90 

DET0318  RDase 

0.62 

DET0145  [Fe]  hydrogenase,  HymA 

0.89 

DET1560  sensory  box  sensor  histidine  kinase 

0.61 

DET0111  [Ni/Fe]  hydrogenase 

0.89 

DET0178  sensor  histidine  kinase 

0.60 

DET0187  formate  dehydrogenase 

0.89 

DET1561  DNA  binding  regulator,  LuxR  family 

0.55 

DET0186  formate  dehydrogenase 

0.88 

DET0316  DNA  binding  response  regulator 

0.51 

DET0343  ftsZ-1  cell  division  protein 

0.87 

DET0305  sensory  box  sensor  histidine  kinase 

0.48 

DET0996  rpmG  ribosomal  protein  L33 

0.87 

DET1540  DNA  binding  response  regulator 

0.40 

DET1225  hypothetical  protein 

0.87 

DET0315  sensory  box  sensor  histidine  kinase 

0.33 

DET0997  tuf  translation  elongation  factor  Tu 

0.87 

DET0322  hypothetical  protein 

0.33 

DET0994  nusG  transcription  antitermination 

0.87 

DET1558  RDase  anchoring  protein 

0.29 

DET1417  purH  cyclohydrolase 

0.86 

DET1164  iron  sulfur  cluster  binding  protein 

0.25 

DET0001  dnaA  replication  initiator 

0.86 

1  Table  4.  Correlation  scores  of  tceA  to  RDases  and  response  regulators  adjacent  to  RDases  and 

2  correlation  scores  of  tceA  to  all  transcripts  on  DET  s  genome  represented  on  the  microarray. 
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1  Figure  5.  RDase  distribution  across  different  electron  acceptors  with  the  (a)  natural  value  ratio 

2  compared  to  a  3-day  starved  culture  for  the  17  RDase  and  the  rpoB  transcript  is  plotted  for  the 

3  PCE,  TCE,  DCE,  DCP,  and  No  Acceptor  fed  conditions.  Error  bars  represent  standard  error 

4  across  biological  replicates,  (b)  The  absolute  processed  intensity  value  distribution  of  the  RDases 

5  for  the  No  Acceptor,  DCP,  DCE,  TCE,  and  PCE  cultures  respectively  moving  inwards  from  the 

6  outer  ring  of  the  pie  chart. 


Analyte  Concentration  (|rmoles/L) 


Transcript  ID  and  Description 

Regulation 

Ratio 

Transcript  ID  and  Description 

Regulation 

Ratio 

Transcript  ID  and  Description 

Regulation 

Ratio 

DET0014  hypothetical  protein 

UP 

2.7 

DET0754  hypothetical  protein 

UP 

2.7 

DET  1 146  hypothetical  protein 

DOWN 

0.36 

DET0035  gmk  guanylate  kinase 

DOWN 

0.42 

DET0755  hypothetical  protein 

UP 

3.8 

DET1 173  fwdE  family  protein 

UP 

2.6 

DET0036  hypothetical  protein 

DOWN 

0.44 

DET086 1  hydrogenase,  EchB  subunit 

DOWN 

0.37 

DET  1186  hypothetical  protein 

UP 

2.9 

DETO 1 08  endonuclease/exonuclease/phosphatase 

UP 

2.2 

DET0862  hydrogenase,  group  4,  EchC  subunit 

DOWN 

0.35 

DET  1206  metM  homoserine  dehydrogenase 

DOWN 

0.34 

DET0137  mgs  A  methylglyoxal  synthase 

DOWN 

0.31 

DET0864  hydrogenase  subunit  HymB 

DOWN 

0.37 

DET  1224  cobA  3  cob(I)alamin  adenosyltransferase 

DOWN 

0.43 

DETO  1 80  reductive  dehalogenase 

UP 

5.5 

DET0908  arsenical  pump  membrane  protein 

UP 

7.8 

DET  1278  nusB  N  utilization  substance  protein  B 

DOWN 

0.25 

DET0181  reductive  dehalogenase  anchoring  protein 

UP 

5.1 

DET0909  hypothetical  protein 

UP 

2.8 

DET  1284  hypothetical  protein 

DOWN 

0.38 

DET0297  hypothetical  protein 

UP 

4.9 

DET0945  ABC  transporter,  permease  protein 

DOWN 

0.43 

DET1285  serine  protease,  DegP/FltrA  family 

UP 

2.9 

DET0298  hypothetical  protein 

UP 

4.5 

DET0964  red  single  stranded  DNA  exonuclease  RecJ 

DOWN 

0.38 

DET  1286  serine  protease,  DegP/FltrA  family 

UP 

3.9 

DET0322  hypothetical  protein 

DOWN 

0.35 

DET1021  hypothetical  protein 

DOWN 

0.20 

DET  13 00  universal  stress  protein  family 

DOWN 

0.27 

DET0368  proS  prolyl  tRNA  synthetase 

UP 

2.7 

DET 1031  hypothetical  protein 

DOWN 

0.40 

DET1377  hypothetical  protein 

UP 

3.0 

DET0420  hypothetical  protein 

DOWN 

0.25 

DET1035  tryptophan  synthase  subunit  beta 

UP 

3.2 

DET  1379  auxin  responsive  GH3  protein  homolog 

UP 

3.7 

DET0423  hypothetical  protein 

UP 

4.2 

DET  103  8  ileS  isoleucyl  tRNA  synthetase 

UP 

2.5 

DET1419  transcriptional  regulator,  AbrB  family 

DOWN 

0.29 

DET0485  rplX  ribosomal  protein  L24 

DOWN 

0.43 

DET  1051  hypothetical  protein 

UP 

3.8 

DET  1 448  hypothetical  protein 

UP 

2.8 

DET0558  atpB  ATP  synthase  F0,  A  subunit 

DOWN 

0.21 

DET1057  hypothetical  protein 

UP 

2.3 

DET151 1  hypothetical  protein 

DOWN 

0.14 

DET0560  atpF  ATP  synthase  F0,  B  subunit 

DOWN 

0.33 

DET  1063  DNA  binding  regulator,  LuxR  family 

UP 

2.7 

DET  1 544  reductive  dehalogenase  anchoring  protein 

UP 

3.5 

DET0561  atpFl  ATP  synthase  FI,  delta  subunit 

DOWN 

0.33 

DET 1068  recombinase,  phage  integrase  family 

DOWN 

0.30 

DET  1 545  reductive  dehalogenase 

UP 

3.3 

DET0562  atpA  ATP  synthase  subunit  A 

DOWN 

0.41 

DET  1092  DNA  methylase 

DOWN 

0.41 

DET1569  hypothetical  protein 

UP 

3.0 

DET0673  hypothetical  protein 

UP 

3.6 

DET  1105  hypothetical  protein 

DOWN 

0.27 

DET  1624  hypothetical  protein 

UP 

2.6 

DET0721  hypothetical  protein 

DOWN 

0.36 

DET  1 1 06  hypothetical  protein 

DOWN 

0.33 

DET0736  oxidoreductase,  dehydrogenase/reductase 

UP 

2.7 

DET  1107  hypothetical  protein 

DOWN 

0.29 

Culture  Normalized  tceA  Processed  Flourescent 
Intensity  (Pfu/mL  Culture/%  DET16S) 
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